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For macroscopic quantum systems, we study what are measured when equilibrium fluctuations of 
macrovariables are measured in an ideal way that mimics classical ideal measurements as closely as 
possible. We find that the symmetrized time correlation (symTC) is always obtained for such mea¬ 
surements. As an important consequence, we show that the fluctuation-dissipation theorem (FDT) 
is partially violated as a relation between observed quantities in macroscopic quantum systems even 
if measurements are made in such an ideal way. 
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When temporal fluctuations of observables are mea¬ 
sured in quantum systems, disturbances by the measure¬ 
ment play crucial roles w- Suppose that one mea¬ 
sures fluctuation of A by measuring A at time t = 0 
and subsequently at time t > 0. Then the disturbance 
by the first measurement affects the outcome of the sec¬ 
ond one. Consequently, fluctuation of A, expressed by its 
time correlation, depends strongly on the way of measur¬ 
ing it. For example, the normal-order time correlation 
and antinormal-order time correlation are obtained, re¬ 
spectively, when a photon field is measured with a pho¬ 
todiode and a quantum counter m • This means, for 
example, that the former cannot measure the zero-point 
fluctuation [3]. 

Such strong dependence on measuring apparatuses 
should also be present in measurements of equilibrium 
fluctuations in macroscopic quantum systems. This may 
affect fundamental relations such as the FDT, which is a 
universal relation between response functions (that char¬ 
acterize responses to driving forces) to equilibrium fluc¬ 
tuations (that are expressed by time correlations) [MU- 
However, the disturbances by measurements were com¬ 
pletely disregarded when deriving the FDT [8UI2]. It is 
obvious that violent measurements of time correlations 
will lead to apparent violation of the FDT because of 
strong disturbances. Hence, one would expect that the 
FDT would hold only when certain ‘ideal’ measurements 
are made for time correlations. For classical systems, 
ideal measurements are those which do not disturb the 
system at all. For quantum systems, however, such mea¬ 
surements are generally impossible. Hence, ‘ideal’ mea¬ 
surements are those which mimic classical ideal measure¬ 
ments as closely as possible. It is expected that the FDT 
should hold when such ‘ideal’ measurements are made. 
Is this expectation true? 

A similar problem exists in Onsager’s regression hy¬ 
pothesis [HI Q3] . While Refs. [T3HIB] claimed its consis¬ 
tency with the Kubo formula m, other works claimed 
inconsistency [T7WTT)1 . To settle this controversy in an op¬ 
erational manner, disturbances by measurements should 
be considered, which were disregarded in these works. 


In this paper, we study, using modern theory of quan¬ 
tum measurement [DEI, what is observed when ‘ideal’ 
measurements are made on equilibrium fluctuations. We 
find the universal answer that the symTC is always ob¬ 
served. As a striking consequence, the FDT and the 
regression hypothesis are violated as relations between 
observed quantities. Furthermore, we show that in the 
post-measurement states, unlike in the Gibbs states, ex¬ 
pectation values of macroscopic observables evolve with 
time. Such states should be realized in experiments on 
equilibrium fluctuations. 

Assumptions on equilibrium state .— We consider a d- 
dimensional macroscopic quantum system (d > 1), whose 
size is characterized by N (such as the number of spins). 
The pre-measurement state, represented by the canonical 
Gibbs state p eq , is its equilibrium state at finite temper¬ 
ature T (= 1//3). 

We assume that the correlation between any local 
observables at two points r and r' decays faster than 
l/|r — r'\ d+t (e > 0) with increasing |r — r'\. This 
assumption is believed to hold generally except in spe¬ 
cial regions, such as critical points, in the thermody¬ 
namic configuration space. Under this assumption, the 
static equilibrium fluctuation, dL4 eq = (((AA) 2 ) eq ) 1 / 2 , of 
any additive observable A is 0{sfN) [2UH221 . Hereafter 
( • ) e q = Tr[p e q • ] and A denotes the shift from the equi¬ 
librium value, i.e., A A = A — (A) eq . Furthermore, with 
some reasonable additional assumptions (see Ref. [23]), 
the ‘quantum central limit theorem’ [2311271 holds, which 
enables us to derive the following universal results. Its 
statement, discussion on related theorems [28H30] , and 
details of the calculations are described in Ref. [23]. 

Quasiclassical measurement. — As discussed above, 
certain ‘ideal’ measurements, which mimic classical ones, 
should be made to get time correlations correctly. We 
call such measurements quasiclassical. To be definite, we 
define them as minimally-disturbing, homogeneous, and 
unbiased quantum measurements with moderate magni¬ 
tudes of errors, as follows. 

To observe fluctuation of an additive observable A, its 
measurement error <JA err should be smaller than <5A cq . 
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On the other hand, <T4 err should not be too small be¬ 
cause otherwise disturbances by the measurement would 
be too large so that the measurement would be much dif¬ 
ferent from classical ideal ones. Therefore, we consider 
measurements in which SA OTV = e5A eqi where e is inde¬ 
pendent of N . Although we are interested in the case 
where e < 1, the following results hold also for larger e 
(which occurs, e.g., when interaction with the measuring 
apparatus is weak my Considering <U4 oq = 0(y/N) as 
mentioned above, we require more concretely (i) the mea¬ 
surement operator (that gives the state just after mea¬ 
surement 00 ) scales as \fN with increasing N. A qua- 
siclassical measurement should also possess the following 
reasonable properties: (ii) Minimally disturbing, i.e, the 
only disturbance of the system is the necessary backac- 
tion determined by the probability operator [5]. (iii) Ho¬ 
mogeneous, i.e., the measurement operator is a function 
of A — A, , where A . is the outcome of the measurement. 
This implies, e.g., a reasonable condition that 8A eTI is in¬ 
dependent of A,, (iv) Unbiased, i.e, A, = (A) eq , where 
7TT denotes average over many runs of experiments. From 
(i)-(iv), the measurement operator should take the form 
f(a — a .) in terms of scaled quantities a = A/y/N and 
a, = A,/y/N. Here, f(x ) is a real function independent 
of TV such that f \f(x)\ 2 dx = 1, f x\f(x)\ 2 dx = 0. Hence, 
the probability density for getting a. is 

P(«.) = <{/(«-a.)} 2 ) eq , (1) 

and the post-measurement state for the outcome a, is 

m 


/3(a.) = /(a - a,)p eq f{a - a,)/p{a,). (2) 


It seems reasonable to require also (v) f(x) behaves 
well enough, e.g., it vanishes quickly with increasing |x| 
(such as the Schwartz functions [551155] : see Ref. [25] 
for details). Although (ii) requires also that / should 
be nonnegative up to an irrelevant phase factor [5], all 
the following results hold without this additional condi¬ 
tion. Then, the measurement error Sa err (= 5A elI /y/N) 
is given by ( Sa erl ) 2 = J x 2 \f(x)\ 2 dx = 0(1), in consis¬ 
tency with 8A eiI = £SA eq . A typical / is gaussian [57], 
f(x) = (27tui 2 ) _1//4 exp(— x 2 /4w 2 ) with w = 0(1) > 0, 
for which Sa elI = w. 

Which time correlation is measured ?— Although all 
the following results hold in the thermodynamic limit, 
we do not write limjv^oo explicitly. Suppose that the 
macroscopic system was in a Gibbs state p eq for t < 0, 
and that an additive observable A is measured at t = 0. 
The outcome a. (= A,/yfN ) distributes according to 
Eq. Q , for which we have 


p(«.) = J 


/' l/( * )|2 c-n 

(x + Aa.) 2 

/ (27T(5a 2 q ) 1 / 2 P 

2Sa eq . 


where 6a 2 q = ( SA eq ) 2 /N and A a, = a, — (a) e 
means that p(a % ) is a convolution of the distribution in 


dx, 

(3) 

') eq ■ 

This 


p eq and the shape |/(:r)| 2 of the measurement operator. 
We also see that the variance scales as Sal — (Aa.) 2 = 
0(1), in consistency with <5a 2 q = 0(1) and <5a 2 rr = 0(1). 

The post-measurement state p(a.) deviates from p eq , 
hence the expectation values, denoted by (•)„., of ob¬ 
servables evolve with time. We investigate the time evo¬ 
lution of an additive observable B (= y/~Nb) for t > 0. 
For (A b(t)) a . = ( b(t)) a . - (b) e q, we find [3F, 5? 

(A b(t)) a . = -0(f)(i{Aa, A6(t)}) eq (lnp) / , (4) 

where 0(f) is the step function, {\{X, F(f)}) eq = 
{\{XY{t) + Y{t)X)) e q the symTC, b(t) = e iflt l h be~ iflt / h 
the Heisenberg operator, and (lnp)' = din p/da,. Since 
p' vanishes at A a, = 0 from Eq. (|3| , the rhs of Eq. 0 is 
linear in Aa« for small |Ao,|. 

Averaging Eq. 0 over a,, we find [3B] 

(A b{t))a, = 0, i.e., (S(t))a. = (b)e q (5) 

for all additive observable B = y/Nb at all t > 0. That 
is, the measurement does not cause any systematic dis¬ 
turbance on b 1 in consistency with our requirements on 
quasiclassical measurements. 

With increasing t, ( b(t)) a . relaxes to (b) eq if the system 
possesses ‘mixing property’ in the sense that 

lim (|{Aa, A6(f)}) eq = 0, (6) 

£—>•00 

where lirri 7 v->.oo is taken before lim^oo. After the relax- 
ation, one cannot distinguish p(a,) from p eq by measur¬ 
ing any additive observables. Equation 0 shows that 
the relaxation process after the measurement is governed 
by the symTC. To see it clearly, we calculate the corre¬ 
lation between a. and (b(t)) am . Multiplying Eq. 0 with 
Aa., and averaging over a,, we find [3B] 

Aa.(Ab(t))a. = 0(t)(|{Aa, A6(t)}) eq . (7) 

This is a universal result independent of choice of f(x), 
and tells us the operational meaning: When one measures 
a time correlation quasiclassically, what he observes is the 
symTC rather than many other time correlations which 
reduce to the same classical correlation as h —> 0. 

This result might look contradictory to some experi¬ 
ments dm mm- However, those measurements are 
not quasiclassical because they destruct the states by ab¬ 
sorbing quanta. If, e.g., heterodyning techniques [331135] 
or quantum non-demolition photodetectors j46fi50| are 
used, the symTC will be observed. 

Violation of FDT.- The above finding has a great im¬ 
pact on nonequilibrium physics. The Kubo formula m 
gives the response function, which describes response of 
a quantum system to an external force, by the canonical 
time correlation, denoted by (■; • ) eq , of certain additive 
observables as [53] 

®ba(t) = 0(t)/3(Aa; Ab(t)) eq . 


( 8 ) 
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Here, the step function 0(t) comes from the causality. 
This formula has been regarded as the FDT for quantum 
systems [T0HT2] . In its derivation, however, disturbances 
by quantum measurements were neglected, although they 
should be considered seriously m■ Therefore, we here 
take Eq. Q just as a recipe to obtain the response func¬ 
tion, as discussed in Ref. [23j , while measured fluctuation 
may possibly be described by a different expression [52] . 
Then a question arises; does the FDT hold as relations 
between measured quantities in quantum systems? 

To answer this fundamental question, we calculate the 
Fourier transform of &b a (t) : which is the admittance de¬ 
noted by \b a (w), and that of the measured time corre¬ 
lation Eq. (J7|, denoted by Sb a (uj) (whereas the Fourier 
transform of the symTC without 0(t) is denoted by 
Sba{u))- We express the results in terms of the sym¬ 
metric and antisymmetric parts, denoted by + and —, 
respectively, such as xta = ( Xba ± Xab)/2. For the real 
parts, we have [53] 


Re X&~ 0 M = £ R 

Re XtTaM = /3Re5' b "(w) 


+P 





1 

WO. 


iSbaiu}') 


cLo' 

2 tT‘ 


(9) 

( 10 ) 


Here, V denotes the principal value, and Ip{ui) = 
(/3hoj/2) coth(/?fiw/2). [For the imaginary parts, replace 

Re Xfia; ReS ba and iSba with Im X'L Im S ba and re¬ 
spectively.] Since Ipipj) ^ 1 at finite lo , we find that the 
FDT is violated at finite uj, when comparing the observed 
equilibrium fluctuation and the observed admittance even 
if the measurement is quasiclassical (i.e., even if it mimics 
classical ideal measurement.) 

One might expect that the FDT would recover in the 
‘classical regime’ where Huj <C k^T, because Ip(vj) -» 1 . 
We examine this expectation by studying the case of w = 
0 (for which y' 6 a ( 0 ) = Re Xi, a ( d ) because &ba(t) is real). 
From Eqs. 0 and ( [To] ) , we find that even at ui = 0 the 
FDT is recovered only for the symmetric part. For the 
antisymmetric part, the causality in Eq. ([ 8 ]) (represented 
by 0 (f)) convolutes different frequencies as 


Xfeq(O) 

P 


-Sbai 0 ) = 



V 

UJ 



1 

J /3M. 


( n ) 


Since (w) is a pure-imaginary odd function [ 12 ], the 
rhs of Eq. © does not vanish in general. Therefore, in 
general, the FDT is violated even in the ‘classical regime’ 
for the antisymmetric part even if the measurement is 
quasiclassical. Note that there are two ways to reach 
the ’classical regime’ flhco -C 1. One is to take h —> 
0 , where the system becomes classical and the violation 
disappears. The other is to take w —>■ 0 while keeping 
h constant, where the violation occurs. Therefore, the 
violation is a genuine quantum effect. 

Experiments .— As an example, we compare the con¬ 
ductivity tensor [23], a, xv {uj) = 8 / 0 °° Q„\Jn(t)) eq e iut dt, 


with the Fourier transform of the measured equilibrium 
fluctuation, S^{uj) = J^° (^{j^Jy.{t)}) eq e lu,t dt, where > 
denotes the v component (v = x , y, z) of the total current 
divided by \[N. 

Koch et al. [2] measured the diagonal (y = v) ele¬ 
ments of the equilibrium fluctuation of a circuit by using 
the heterodyning technique [45] ■ Their measurement is 
closer to quasiclassical than those of Refs. [38H43] be¬ 
cause it does not destruct states by absorbing quanta, 
as discussed above. Hence, it can be regarded as a pio¬ 
neering work about the equilibrium fluctuation obtained 
by quasiclassical measurements. Its results are consistent 
with Eq. ([£]). 

For the off-diagonal (y 7 ^ u) elements, no experiments 
have been reported to the authors’ knowledge. We pro¬ 
pose to examine the w = 0 case, Eq. (11), by studying 
<T xy { 0 ) in the presence of a magnetic field parallel to the 
z axis. When the system is invariant under rotation by 
7t/2 about the z axis, the obvious symmetry a xy = —a yx 
and Eq. (11) yield 


P 




1 


iS"» —. (12) 


Therefore, even at cu = 0 the FDT is violated for the 
off-diagonal (Hall) conductivity. This violation will be 
confirmed by measuring independently cr xy (0) and the 
correlation jy(j x (t))j ym , which gives S xy (ui) according 
to Eq. ([7]). One can also obtain the rhs of Eq. ]l2[ by 
measuring and which together give S^u) 

according to Eq. © below. 

Regression hypothesis .— References IMS! pointed 
out that Onsager’s regression hypothesis msui] would 
contradict with Eq. ] 8 |, assuming the symTC for the 
time correlation of the hypothesis. Nakajima showed that 
this contradiction can be removed if a local equilibrium 
state is assumed for the state during fluctuation [15] , and 
thereby derived response to non-mechanical forces [16] . 
These contradictory claims were derived from different 
assumptions, but none of them was justified satisfacto¬ 
rily. Since we have proved that the symTC is always 
measured in quasiclassical experiments, the hypothesis 
cannot be valid as a relation between measured quanti¬ 
ties in quantum systems even if the measurements mimic 
classical ideal measurements. 

Squeezed equilibrium state .— We have shown how 
quantum effects appear in the expectation value, Eqs. 0 
and ([5]), and in the correlation, Eq. (J7|. Quantum effects 
appear more manifestly in the variance of b(t), which is 
calculated as m 


mt)-(b(t)) a .) 2 ) a .-6bi q 


= (k{Aa,Ab(t)}) 2 eq (ln P r- 


(hMm. 


2 

eq 


d 2 \a.q 

dy 2 


,(13) 

v=0 
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where 


q(a.,y) = (f(d - a. + y)f{d -a,- y)) eq 


[ f{x + y)f{x-y) 

(x + Aa.) 2 

1 (27n5a| q ) 1 / 2 ' l> 

2Sa eq . 


dx. (14) 


We find that the relaxation process is governed by 
both the symTC and the commutator time correla¬ 
tion (7^ [a, b(t)]) eq . If, in addition to Eq. (JHJ) , the sys¬ 
tem also has the mixing property in the sense that 
lim^oo limAr_ ) . 00 (^[a, 6(t)])g q = 0, then the variance re¬ 
laxes to 6b 2 q with increasing t. 


We can see the physical meaning of Eq. (13) by let 


ting t —> +0. The symTC implies that p(a,) is squeezed 
along b by measuring a if b correlates with a in p eq . The 
commutator time correlation implies that the squeezing 
is disturbed by the measurement if [a, b] is non-negligible 
in /5 eq . Furthermore, the rhs of Eq. Bl) is 0(1), which is 
the same order as i5o eq . That is, p(af) deviates from p eq 
only within equilibrium fluctuations, hence the system re¬ 
mains nracroscopically in the same equilibrium state. In 
this ‘squeezed equilibrium state,’ macrovariables evolve 
with time as Eq. unlike in the Gibbs state. It rep¬ 
resents the state that is observed in quasiclassical mea¬ 
surements of equilibrium fluctuations. 

Multi-time measurements. — To measure the lhs’s of 
Eqs. Q and ([ 7 ]), one must perform measurements twice 
in each run of the experiment, as described in Ref. [23] - 
When one performs more measurements subsequently in 
each run, what he gets is the following. 

Suppose that K+l additive operators A 0 , A 1 , - ■ ■ , A K 
are measured subsequently at t = to, ti, ■ • • , tx, re¬ 
spectively, and the outcomes are A^,A\,--- ,A, (= 

\/Na ®, y/~Na \, • • • ,y/Najf). Here, K = 0(1) and 0 = 
to < t\ < ■ ■ ■ < tx- Let the measurement operator for 
Ai = y/Na? be /j(a J — aj), where fj(x) satisfies the 
aforementioned conditions for f{x). The state at t = tj 
is obtained by applying fj(d j —ai) ■ ■ ■e~’ l ^ I( d 2 ~ tl ' > / h 
a\)e~ lHtl / h fo(a°— a 0 ,) and its conjugate to p eq from the 
left and right, respectively. We find that 


A ai = 0, i.e., aj = (a J ') eq for all j. (15) 

For the correlations for 0 < j < k < I \, we get 

Aa 3 .Aa k = (l{Ad J (t j ),Ad k (t k )}) eq + Sj,k5a J J r 

3 -1 

+ < Ti ^ tty ) > &l Ob)] )eq 1 {ji [a 1 (ti) , d k (4)]) eq ., (:ifi) 

;=0 

where = J x 2 \fj\ 2 dx and Fj = 4 f {f'(x)} 2 dx 137] . 
The first term in the rhs of Eq. © corresponds to the 
correlation in Eq. ([ 7 ]). The second term comes from 
the measurement error of each measurement, which is 
absent in Eq. 0 because it corresponds to the case 
Sj t k = <5q,i = 0- The last term represents disturbances 


by the measurements that are performed before the j-tii 
measurement, hence the summation is over m such that 
m < j. Because of this term, which depends on the ex¬ 
perimental setup represented by {dj, fj(x),tj}j, the cor¬ 
relation deviates from the symTC unlike the case of twice 
measurements, Eq. ([7]). Consequently, the FDT is vio¬ 
lated more drastically in this protocol of experiment , e.g., 
even at uj = 0 both the symmetric and antisymmetric 
parts can violate the FDT. 

For the special case where j = 0 (and k > 1), however, 
the disturbance term is absent: 

Aa°Aa^ = (L{Aa°, Ad fc (t fe )}) eq for t k > 0. (17) 

This coincides with Eq. ([T]), although other measure¬ 
ments may be performed for 0 < t < t k - It is also univer¬ 
sal, independent of choice of f{x), as is Eq. Q. Hence, 
in this case, the FDT is violated in the same manner as 
in the previous protocol. 

In summary, we have studied what is measured when 
the equilibrium fluctuation is measured in an ideal way 
that mimics classical ideal measurements as closely as 
possible, i,e., such that disturbances are as small as pos¬ 
sible under the condition that the equilibrium fluctuation 
can be measured accurately enough. We call such mea¬ 
surements quasiclassical. It is found that the symTC 
is obtained quite generally [Eqs. ([T]) and ©]• From 
this finding and the causality, we have shown that the 
FDT is violated as a relation between observed quan¬ 
tities even if measurements are quasiclassical [Eqs. ®>- 
jnj)]. It is violated for antisymmetric parts of response 
functions at all frequencies and for symmetric parts at 
high frequencies hu> > fc^T. Onsager’s regression hy¬ 
pothesis is also violated. These striking results are gen¬ 
uine quantum effects. Disturbances by measurements 
appear more strongly in the variances of macrovariables 
[Eq. ( T 3 I)] and in the case of multi-time measurements 
[Eq. (Ifif]. The post-measurement state is shown to be 
a time-evolving ‘squeezed equilibrium state,’ in which 
macrovariables fluctuate and relax [Eqs. <@ and ©] . It 
represents the state realized during quasiclassical mea¬ 
surements of equilibrium fluctuations. 

Finally, we note that our results apply not only to the 
Gibbs states but also to pure states representing equilib¬ 
rium states [53MS7] , 
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We first present generalization of measurement operators in Eq. (2). We then give a brief review 
of the Kubo formula in Sec.II. There, we also discuss irrelevance of disturbances by measurements 
when measuring response functions. In Sec. Ill, we explain the experimental protocol for measuring 
the lhs’s of Eqs. (4) and (7). In Sec. IV, we show a version of the conditions and the statement 
of the quantum central limit theorem (QCLT), which is used for deriving the results of the paper. 
Furthermore, we present other theorems and their proofs to derive the main results of the paper. In 
Sec. V, we present derivations of the main results of the paper, showing assumptions on the system 
and sufficient conditions on f(x). Finally in Sec. VI, we present the results for the case of a gaussian 
/(*)■ 


I. GENERALIZATION OF MEASUREMENT OPERATORS IN EQ. (2) 

We can generalize Eq. (2) as [1] 

L 

p(a.) oc (1 /L) Y / A (« - a.)p eq f x {d - a,) (LI) 

A=1 

with L different functions f 1 , f 2 , ■ ■ ■ , f L that satisfy the conditions for / of the paper. As long as L = 0(1), one can 
obtain the results for L > 1 by averaging the results of the paper over all f x . 

II. BRIEF REVIEW OF KUBO FORMULA 
A. Kubo formula 

Consider a macroscopic quantum system, which was in an equilibrium state p eq = e~^ H /Z in a remote past. 
Suppose that a spatially-varying time-dependent external field F(r,t) is applied, which interacts with the system via 
the interaction Hamiltonian 


H ex t(t) = ~Y F ( r ’ *)7( r )> (II.l) 

r 

where y(r) is a ‘local operator’ that acts on a finite number of sites around site r of the system. Then the state p(t) of 
the system deviates from p eq . Consequently, the expectation value ( a(r)) t = Tr[p(f)&(r)] of a local observable a(r) 
deviates from its equilibrium value (a(r)) eq = Tr[/3 eq a:(r)]. Their difference is the response of d(r) to F(r,t). To the 
linear order in F(r,t), the response and F(r,t) is related by the (linear) response function 4>(r,t) as 

{a{r))t-(a(r)) eq = f dt' Y <K r - r ', 1 ~ ( IL2 ) 

■' — oo r , 

Although the r dependence is often crucial, 1 we shall present a simplified case where F(r,t) is independent of r, 


* fujikura@ASone.c.u-tokyo.ac.jp 
' shmz@ASone.c.u-tokyo.ac.jp (contact author) 

1 Let k and uj be the wavenumber and frequency of F(r,t). In general, when one is interested in the limit of k —> 0 and cj —> 0, the two 
limits should be taken in an appropriate order to get a correct result. This means, for example, that the result for k —> 0 often disagrees 
with the result for k = 0, hence the r dependence cannot be neglected. See Ref. [3] on this point and other issues. 
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i.e., F(r,t ) = F(f). Then, Eqs. (II.1) and (II.2) are simplified to 

^ext(i) = -F(t)C, 

= *(* ~ t')F(t')dt', 


(IL3) 

(II.4) 


where 


c , = E^ r )’ ( IL5 ) 

r 

^ = ( n - 6 ) 

r 

are additive observables, N is the number of sites of the system, and 

^{t)= ^^2^2(j>{r, r ',t) (II.7) 

r r' 

is the response function describing the linear response of A/N to F(t). 

Since we consider the linear response, F(t) is assumed to be small and independent of N. Hence, ( A) t /N is also 
small and independent of N (for sufficiently large N). Therefore, $(f) = 0(1) according to Eq. (II.4). 

Assuming the unitary time evolution (see Sec. IIC) by the total Hamiltonian FI + ff ext (t), Kubo [2] showed that 


$(f) = 0(f) lim ^([c,a(f)]) eq , 
n^-oo in 

= 0(f) lim /?(&; a(f)) eq , 


(11.8) 

(11.9) 


which is called the Kubo formula. Here, c = C/y/N, a(t) is the Heisenberg operator (e lHt / h ae lHt / h ) of a = A/y/N, 
and b = B/y/N , where 




t=0 


= '■ 


(n.io) 


If H includes only short-range interactions, B becomes an additive observable. Furthermore, ( ; ) eq denotes the 
canonical correlation, which is defined for two operators X and Y by 


<^}eq 



{e x6 X^e~ xft Y) e(l dX. 


( 11 . 11 ) 


Note that 


(X;i) eq = (i;X) eq = (A) eq . (11.12) 

Since 

(B)e q= |(0(t))eq 

we can rewrite Eq. (II.9), using Eq. (11.12), as 

®ab(t) = 0(f) lim f3(Ab; Aa(f)) eq , (H-14) 

N—>oo 

where Ab := b — (b) eq (= b) and Aa(f) := a(f) — (a) eq , and the subscript ‘ ab ’ of 4> a f,(t) denotes the dependence of the 
response function on b and a. This is Eq. (8) of the text (where the roles of a and b are interchanged). 


= 0 , 


t—0 


(11.13) 
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B. Admittance 


The admittance Xabflj) is the Fourier transform of <f> a f,(t). Its symmetric and antisymmetric parts are defined as in 
the text by 


,±, ^ Xab(u) ±Xba(u) 

) ' o 


(11.15) 


Here, Xbaflj) is the admittance for the case where the roles of a and b are interchanged. One is often interested in 
the case where F(t) and A are some components of vectors. In such a case, Ad and A b in Eq. (11.14) become the 
corresponding components of vectors, and the admittance Xab(bj) becomes a tensor, and Xab( u ) symmetric and 
antisymmetric tensor components. 

For example, one can take F(t) = E y (t) (the y component of an electric field E(t)) and A = J x (the x component 
of the total current J). Then, Xab{cj) is the Hall conductivity cr xy (u>). In particular, when the system is invariant 
under rotation by 7r/2 about the z axis, we also have X~ b (w) = < T xy {u >) because of the obvious symmetry cr xy = —a yx . 


C. Irrelevance of perturbations by measurements when measuring response functions 

The FDT claims that 


response function = j3 x equilibrium fluctuation. (11.16) 

Let us consider how accurately the response function and the fluctuation should be measured to confirm the FDT. 

To measure the response function, one applies F(t) to the system, which induces a change of A of O(N). By mea¬ 
suring this change, one can determine the response function $(t) according to Eq. (II.4). Therefore, for measuring the 
response function it is sufficient to measure (A) t with measurement error of O(N). More concretely, the measurement 
error for A should be e(A) eq , where 0 < e <C 1. 

On the other hand, to measure the equilibrium fluctuation, one does not apply F(t) and watch fluctuations of A and 
B. Since their fluctuations are of O(VN), one needs to measure A and B with measurement error of 0(y/N). More 
concretely, the measurement error for A should be £(L4 eq , as stated in the text, where 0 < e <C 1. This measurement 
error is 0(l/y/N) times smaller than that for the measurement of the response function. This means that disturbances 
by measurements affect the fluctuation measurements much more strongly than the measurements of response functions, 
because, in general, the (minimum) disturbance by measurement becomes larger as the measurement error is decreased. 

In fact, we have shown in our paper that disturbances of quasiclassical measurements are 0(y/N) for additive 
observables. 2 For sufficiently large N, such disturbances are negligible for measurements of ( A) t /N although they 
play crucial roles in the fluctuation measurements as our paper has shown. Therefore, disturbances by quasiclassical 
measurements can be neglected in measurements of response functions. For this reason, we have assumed in our 
paper that the Kubo formula is correct as a recipe to obtain the response function (while measured fluctuation 
may possibly be described by a different expression), although Kubo [2] calculated the response function without 
considering disturbances by measurements. 


D. Retarded Green function 

Among two forms (II.8) and (11.14) of the Kubo formula, formula (11.14) is preferred in statistical mechanics 
because it resembles the formula for thermodynamic susceptibilities (i.e., response of equilibrium values to equilibrium 
parameters such as a static magnetic held) for quantum systems. This similarity is convenient, e.g., when a necessary 
condition for the validity of the Kubo formula is discussed [2, 3]. 

On the other hand, formula (II.8) is generally preferred in condensed-matter physics because the rhs is nothing but 
the retarded Green function, which can be calculated using the diagrammatic expansion. In practical calculations, 
formula (II.8) is often transformed into other forms using, e.g., the charge conservation. For example, the longitudinal 
conductivity is usually calculated after transforming formula (II.8) into the current-current retarded Green function. 


2 


Although we have shown this fact for equilibrium states, we expect that the order of magnitude does not change in the linear-response 
regime. 
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III. PROTOCOL FOR MEASURING THE LHS’S OF EQS. (4) AND (7) 

The protocol for measuring the lhs’s of Eqs. (4) and (7) is as follows. 

One must perform many runs of experiments, which we label r = 1,2, • • •. In each experiment, 

1) prepare the system in the same equilibrium state, 

2) measure a at time t r to get the outcome a, r , 

3) measure b at time t r +1 to get the outcome b, r . 

By repeating l)-3) many times, one gets a set of data {a #r , b, r } r . Take arbitrarily a value a,, and select all data that 
satisfy \a, r — «,| < 1. By averaging b, r over the selected data, one obtains ( b(t)) a . ■ By subtracting ( b) eq , a way of 
getting which is obvious, one gets the lhs of Eq. (4). By multiplying this with Aa«, and by averaging over a», one 
gets the lhs of Eq. (7). 


IV. THEOREMS 

We write N for the size, such as the number of spins, of the quantum system. The thermodynamic limit (TDL) is 
briefly denoted by lim throughout this Supplemental Material (and the paper). The justification of interchanging 

N—> oo 

the order of the TDL and differentiation is one of the keys to the rigorous derivations of our results. 

We let 

4 = £&<(»■) (iv. 1) 

r 

be an additive observable, where &i(r) is a local observable acting around r, and r runs all over the volume of the 
system. The index i = 1,2,..., M labels such additive observables of interest, where M = 0(1). 


A. Conditions for QCLT 


The conditions for the QCLT are summarized as follows. 

(Cl) The following quantities are 0(1) (i.e., they have IV-independent upper bounds): 

|[i i ,d i (r)]|| , | [Ai, [Aj,a k (r)] I 


(IV.2) 


(C2) The following quantity decays faster than 1 /|r — r'\ d+e (e > 0 is some constant) with increasing the distance 


r — r . 


(-{Ad,:(r), Ad^r'Deq 


(C3) For any simply connected subsystem A, 


(IV.3) 

(IV.4) 


T. |(Ad i (ri)Ad i (r2)Adfc(r 3 )Ad/(r4))| = 0 (|A| 2 ). 

r\ ,T2,T’3 ) '^4GA 

This implies, for example, that the forth moment of any additive operator is 0(N 2 ). 

(C4) For some positive constant s, 

N 1 / 2 c N (N 1 / 2d ~ £ , jv 1 / 2d+£ ; f0 (N -4 oo), (IV.5) 

where 

c N (r, V';£) 

= supj (e^Af-^'/Wv ... e iZ M A A + A '/s/N}^ _ ^a^/Vn ... e H M A% /vW... e . 

d( A, A) > r,min{|A|, |A'|} < N'\. (IV.6) 
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Here, A and A' are arbitrary simply-connected subsystems, 


and Aj denotes the restriction of Aj to A, i.e., 


A f = MM 

A 


(IV.7) 


B. Quantum Central Limit Theorems 


Theorem 1 (QCLT for the characteristic function). If the conditions of Sec. IVA are satisfied, then 

Jim (expIi-^Aij • • • exp[iAA M ]) eq = exp]-* V fjf k Sj k \, 

Af->oo VN \/N 2 

3->K 

where 

[j im -^(AA,Ai fc ) eq (j<k) 

S jk = )n^ooN 

iiv 1 Tooiv (AAfeAAi)eq {j>k) 

= L 1 ^^i A M A M» eq + ^t A M Ai *»- O' < *0 

[nZo ^\ {AAj ' Aife })eq + ^ ^ Ai ^>«i 0 > *0 


(IV.8) 


(IV.9) 


(IV.10) 


This theorem implies that the normalized fluctuation operators A Aj/y/N behave as c-numbers, i.e. as bosons, in 
the TDL. Indeed, one can consider that, in the TDL, the fluctuation operators A Aj/y/N form the CCR (Canonical 
Commutation Relation) algebra and the state is represented by a gaussian state for the algebra. This theorem is an 
extension to quantum mechanics of the (classical) central limit theorem, which states that the probabilistic distribution 
of an additive quantity tends to obey the gaussian distribution. For the detail and proof of the above theorem, see 

Ml- 

Actually, the conditions of Sec. IV A are a bit stronger than required for the QCLT of this form. The conditions, 
though, are required for Theorem 3. 

Furthermore, we can prove the following QCLTs for the first and second derivatives of the ‘characteristic function,’ 
which is a quantum analog of the characteristic function. These theorems will be used in Sec. V. 

Theorem 2 (QCLT for the first derivative). If the conditions of Sec. IVA are satisfied, then 


lim ^-(exp^-^LAAi] 
N^oo 8^^ ^ y/N 


' ex P[*~ 7 =AA M ]) eq = Jim (exp^^Aij 
VN dfj n^oo vCV 


•exp[i-^Ai M ]) eq 


d_ 

dfj 


exp 


r. J ] fkfl s kl 


k,l 


(IV.ll) 
(IV.12) 


Theorem 3 (QCLT for the second derivative). If the conditions of Sec. IVA are satisfied, then 


jfe, I"* (exp[i 7y' 4 ' 4 ’ 1 ’'' eXP[i 7=F aiMl) ” = ^ ^(txpfi^jAi.1 ■ - -exp[,i| Ai„]). q 



1 

2 


TMish 

k,l 


(IV.13) 
(IV.14) 


To prove these theorems, we prepare some lemmas. Hereafter, we call an operator [/(£) a unitary operator generated 
by fluctuation if 


f/(£) = g*fiAAi /VN _ _ _ e i£ L AA L /VN 

where A\,... ,Al are additive operators, and L is some constant of 0(1). We write |£| = JV \. 


(IV.15) 
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Lemma 1. Assume that the system satisfies condition (Cl). Let [/(£) be a unitary operator generated by fluctuation, 
and B an additive operator. Then, there exists a positive constant I\ of 0(1) for which the following inequality holds 
for any 


1 

VN 


[U,B] < AT|€|. 


(IV.16) 


Proof. First, we prove the lemma for the case U (£) = e < A ^/vN_ 


s/N 


jZaA/Vn 


B} 


1 

N 


f A\e itL ' AA/VN 

Jo 


<1% [A,B] 


By condition (Cl) the rhs is 0(1). Next, we consider the general case where U(fi) = e l ^ Al /'' A ■ ■ ■ e’OWvW. 


Vn 


l e 


iZiM/y/N ... p ie. L A L /y/N ^|| < V he^Aj/^N 


.7 = 1 


,B] 


<*X>I [ApB] 


3 = 1 


Again, the rhs is 0(1) by condition (Cl). 


□ 


Lemma 2. Assume that the system satisfies condition (Cl). Let [/(£) be a unitary operator generated by fluctuation, 
and Bi, B 2 be additive operators. Then, there exist positive constants Ki, K 2 of 0(1) for which the following inequality 
holds for any 


N 


[[M r],S 2 ] < 


(IV.17) 


Proof. First, we prove the lemma for the case U(fi) = e *£ Aj4 /v-V 
iiAA /^ ,B 1 ),B 2 ]\ 


1 

N 


1 


< 


N 3 / 2 

1 

TV 3 / 2 


' ie> i(,t-t,')AA/sfN 


dfiie 


[Bi,A], 


ip AA/y/N 


,B 2 ] 


[ d^’i[e i ^' )AA/ '^,B 2 ][B 1 ,A]e ii ' AA/VN 

Jo 


AT3/2 


d£fie i{ ^' )AA/ ^[[B^ A], B 2 }e^' AA/VN 


TV 3 / 2 


dff ie i( ^' )AA/V ^[13 ll A\[e i ^ AA/VN , B 2 

Jo 


< \\(A,B 2 ] I ||[Bi,i]| + E^ jf df | [[£1,4 S 2 ] | + E J*#' ||[Bi,i]|| \e\\\[A,B 2 ] 

c?± 

- 2 N 2 


[A, B 2 ] (B 1 ,A] [[Br ,A],B 2 ] 


The first term of the rhs is 0(1) and the second term is o(l) by condition (Cl). Next, we consider the general case 
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where U (£) = ^CAi/Vn • • • e l irA L /VN . 

i L „ - .ok. . . 

< _L ^ yCA/ViV ... ( -i, l, /v'A'yi, 1,/v ^'. - I/'/A'... e i^ L A L /VN^ y 2 j 

^ i=i 

i L r, . 

<—^ [e^iAi/VS ... gilj-iAj-x/vTv^^j ^ 

^ j= i L 

+ |[[e^^/^,B 1 ], J B 2 ]| 

_l_ [y^Ai/vT^y^'+i^j+i/vTv _ _ _ 

By Lemma 1 and the previous case, Lemma 2 is proved. □ 

Lemma 3 (Bound for the first derivative of the ‘characteristic function’). Assume that the system satisfies conditions 
(Cl) and (C2). Let Ui(£), {7 2 (£') be unitary operators generated by fluctuation, and B an additive operator. Then, 
there exist positive constants Kq and K\ of 0(1) for which the following inequality holds for any £, 

-^(LMOA^Oeql < Ko + Ktfl (IV.18) 

Proof. 

-T |(£M{)AB!7 2 ({')>« 1 | < -L|(AB£>i(0t/ 2 (S')). q | + -L|{[C,(«),AB]C 2 (f))„| 

<^<AB^ + ^||[& l( € ) .AB]|. 

The last line follows from the Cauchy-Schwarz inequality. The first term is bounded by a constant of 0(1) by condition 
(C2). Also, the second term is bounded by |£| times a constant of 0(1) by Lemma 1. □ 

Lemma 4 (Bound for the second derivative of the ’characteristic function’). Assume that the system satisfies condi¬ 
tions (Cl) and (C2). Let Li(£i) ? 17 2 (£ 2 ), tl 3 (£ 3 ) be unitary operators generated by fluctuation, and B 1 , U 2 be additive 
operators. Then, there exist positive constants Kq, K\ and K 2 of 0(1) for which the following inequality holds for 
any £ 1 , £ 2 , 

^|<tf 1 (£i)AB 1 tf 2 (£ 2 )AB 2 tf 3 (£ 3 )>e q | < Ko + Ad(|£i| + |£ 3 |) + tf 2 |€i||&|. (IV. 19) 

Proof. 

^|(^ 1 (£ 1 )AB 1 f/ 2 (£ 2 )A J B 2 t/ 3 (£ 3 ))e q | 

< ^|(AB 1 I7i(£i)f/ 2 (£ 2 )f/3(6)AB 2 ) eq | + ^|([LTi(£i),AB 1 ]172(6)^3(^3)AB 2 )e q | 

+ y|(AB 1 0 1 (£ 1 )17 2 (£ 2 )[A J B 2 ,?73(l3)])e q | + y 1 (^( 6 ), ABflUifaflAfy, t/ 3 (£ 3 )]) eq | 

< ^(ABDlfiABDlf + 1 ||[?7i(£i),AB 1 ]|| (A B^ 2 

+ y(A Bl)\g |[AB 2 , 17 3 (£ 3 )]|| + 1 ||[E7i(£ 1 ),AB 1 ]|| ||[AB 2 , U 3 (&)\| • 

The last line follows from the Cauchy-Schwarz inequality. From condition (C2) and Lemma 1, Eq. (IV.19) holds. □ 

Lemma 5 (Bound for the third derivative of the ‘characteristic function’). Assume that the system satisfies conditions 
(Cl), (C2) and (C3). Let t/i(£i), t/ 2 (£ 2 ), 17 3 (£ 3 ), Uflfi 4 ) be unitary operators generated by fluctuation, and B±, _B 2 , 
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£>3 be additive operators. Then, there exist positive constants Kq, K\, K 2 and K 3 of 0(1) for which the following 
inequality holds for any £ 1 , 6> 6 and 6-' 

jy 3 / 2 (#l(6)A.Bl# 2 (6)A.B 2 [)3(6)A.B3#4(6))e q 

< Ko + K^l + | 6 | + 161) + AVI 6 II 6 I + I 6 II 6 I + I 6 II 6 I) + -MI 6 I + 161)161161- (IV. 20 ) 

Proof. 

7 ^|([> 1 (6)AB 1 ^ 2 (6)AB 2 ^3(6)AB3C>4(6))e q | 

< I < ABi AB2U1 (6 ) #2 (6 ) #3 (6 ) #4 (6 ) AB 3 )eq I + ^|(AB 1 [[/ 1 (6)^2(6), A4]J7 3 (6)^4(6)AB 3 )e q | 

+ ^l([f>i(6),AB 1 ]^ 2 (6)J7 3 (6)C>4(6)A-B 2 AB3)e q | + ^|([f/l(6),AB 1 ]^ 2 (6)[A4^3(6)^4(6)]AB3)e q | 

+ ^|(AB 1 AB 2 ^i(6)^2(6)^3(6)[A-B 3 ,c/4(6)])e q | + | (AS! [Ui (6) t/h (6), ab 2 ] u 3 (6) [ab 3 , u A (6)] ) eq I 

+ ^|(A4[^ 1 (6),AB 1 ]^ 2 (6)%(6)[A4^4(6)])e q | + ^|([[C/l(6),A J B 1 ],AB 2 ]^2(6)^3(6)[AB3 ) f/4(6)]>e q | 

+ ^l<[^l(^l)> A ^l][^(6),AB 2 ]f/3(6)[A J B3^4(6)])e q | 

< ] ^ 7 2(AB 1 AB 2 2 AB 1 )V2(AB|)V2 + JL ||[J7 1 (6)f/ 2 (6),A J B 2 ]|| (AB?)^{A4 2 )^ 2 

+ ^ ||[^i«i),AB 1 ]|| (ABgABlABa)^ 2 + ^ ||[#i(6),AB 1 ]|| ||[AB 2) # 3 (6)# 4 (6)]|| (AB 3 2 )ef 
+ ^(ABASIAS!)V 2 | [AJ 5 3 , C 7 4 ( 6 )] I + ^ || [#i(6)^(6), AS,]|| |[AB 3 , *6(6)]|| (AB?)^ 2 
+ ^ |^i(^)’AS 1 ]| |[A J B 3 ,t> 4 ( 6 )]| (AB 2 )ef + ^72 ||[[^i(6), AS,], AB 2 ]|| |[AB 3 , # 4 (6)]| 

+ ^ 1& (&)’ A ^ l] | | [^(6), ab 2 ] 11 [AB 3 , &(6)] | - 

The last line follows from the Cauchy-Schwarz inequality and the definition of the operator norm. By conditions (C2) 
and (C3) and Lemmas 1 and 2, Eq. (IV.20) is proved. □ 

Lemma 6. Let I be a closed and bounded interval, and {gN(x)}N be a sequence of twice differentiable functions on I 
which converges (pointwise) to g(x). Assume that the sequence of second derivatives {g'ff(x)}isr are uniformly bounded 
on I. Then, {g'jy(x)}N converges uniformly to g'(x) (hence we can interchange the differential operation and the limit 
operation). 

Proof. We use a corollary of the Arzela-Ascoli theorem : If a sequence of functions {h]y(x)}N on a closed and bounded 
interval is uniformly bounded and equicontinuous, and if any subsequences (of the sequence) which converge uniformly 
have the same limit function h(x) independent of the choice of the subsequence, then the original sequence {Iin(x)}n 
also converges uniformly to h(x). Now, by the condition, the sequence of the first derivatives {( ar )uniformly 
bounded and equicontinuous on I. Additionally, if a subsequence {96 (a;)}*, (Ni < N 2 < ■ ■ ■) converges uniformly, 
its limit is g'(x), independent of {-/V*,}*,. Therefore, { 9 ^( 2 :)}^ converges uniformly to g'(x). □ 

Now, we can prove the theorems using these lemmas. 

Proof of Theorems 2 and 3. We fix f and a closed bounded interval I including £. By Lemma 4, the second derivative 
of ( e l ii AA i/vN ... w ith respect to f t is uniformly bounded on the interval. We can then apply Lemma 

6 , hence Theorem 2 holds. Similarly, by Lemma 5, the third derivative of ( e l £i AA i/vN ... e l ^ lAAl ^'^) eq with respect 
to is uniformly bounded on the interval, hence Theorem 3 holds. □ 

C. Some more theorems 


In this subsection, we present some more theorems which will also be used in the next section for deriving our 
results. The theorems will be proved under certain conditions on f(x), which will be sufficient conditions for our 
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results. We also prove that Schwartz functions satisfy all the conditions. As in the main paper, we write 

a = A/VN , (IV.21) 


A a = a- (a) eq . 


(IV.22) 


Note that a function of a self-adjoint operator, such as f(a — a,), is defined in terms of the spectrum decomposition 
of the operator. 

Theorem 4 (‘characteristic function’ of squeezed equilibrium state). If the conditions of Sec. IVA are satisfied for 
Aq, ..., Am,B, and if fo, ■ ■ ■, fM are bounded piecewise-continuous real functions, then 


lirn (/ 0 (a° - a°) • • • f M {a M - a™)e luAb f M (a M - a • • • / 0 (a M - a^)) eq 
/o[£]/o[£'] • • • /m[£]/m[£'] 


N—>oo 

f dtodd'o d^d^' M 


2lT 


2n 

-0 On 


x lim _a ^) e ^A6 e z£^(a M ~ a ™) ... e ^o( a ° _a 2)\ 

N—±oc> G 


(IV.23) 


dtodt'o d^Md(' M 


2i r 


x exp 


2tt 

M 


/oKl/oKI • • • /m[£]/m[£'] 


-i £(& + &)Aa* - I + &)& + &') limJ-^AaAtA^a 1 ^)}) 

k —0 

1 ^ M 1 

9 u2 - I A a fc (4)}) eq 

/c—0 


2 N—>oo 
M 


M 


- ~ &)(&+$) J im (^[a fc (4),a'(t/)])eq - « -£fc) lim <^[o fe (*fc), *])e 

z —' N—>oo Z z —' N—>oo Z 

k<l k —0 


N—>oc 2 

Here, fk[f] is the Fourier transform of fk(x): 

/*[£] = 




V2n 


fk(x). 


(IV.24) 


(IV.25) 


Proof. By Theorem 2 in [7], the limit of expectation value in the lhs of Eq. (IV.23) is determined by the QCLT, and 
the limit can be expressed by the Fourier transform as (IV.24). □ 

Theorem 5 (expectation value in squeezed equilibrium state). Assume that the conditions of Sec. IVA are satisfied 
for Aq, ..., Am, B, and that fo, ■ ■ ■, /m are bounded piecewise-continuous real functions. If 

(/o(a° — a°) • • • /m(o m — a^ I )Ab 2 fM(a M — a■ ■ ■ fo(a M — a^ / )) eq < an upper bound independent of N, (IV.26) 

then 


lim (/ 0 (a u - a") • • • f M (a M - a™)Abf M (a M - a™) ■ ■ ■ f 0 (a M - af )) eq 

iv—>-oo 


d 


= - a.) • • • /m(o - a. )e lu f M {a - a, ) • • • f 0 (a - a . )) eq | u=0 

Proof. By the Cauchy-Schwarz inequality and (IV.26), 

</o(o° - a°) • • • f M [a M - a^)Ab 2 e iui f M (a M - af) • • • / 0 (a M - af )) eq 

< (fo(a° - a°) • • • f M (a M - a™)Ab 2 f M {a M - a") • • • / 0 (a M - a. M )) eq 

< an upper bound independent of N. 

This implies that the second derivative of the ‘characteristic function’ 

d 2 


(IV.27) 


^ 2 </o(o° -<)■■■ f M (a M - a™y ub f M (a M - a ™) • • • / 0 (« M - af )) eq 


(IV.28) 

is uniformly bounded on any bounded interval of u. Hence, by Theorem 4 and Lemma 6, Eq. (IV.27) is proved. □ 
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Theorem 6 (variance in the squeezed equilibrium state). Assume that the conditions of Sec. IVA are satisfied for 
Ao ,..., Am, B, and that fo, ■ ■ ■, are bounded piecewise-continuous real functions. If 

(/o(a° — a®) • • • /m(o M — a^)A6 4 /m(o M — a^ 7 ) • • • fo{a M — a^ f )) eq < an upper bound independent of N, (IV.29) 

then 


lim (/ 0 (a° - a°) ■ ■ ■ fM{a M - a™)Ab 2 f M (a M - a■ ■ • / 0 (a M - a 7/ )) eq 

Af->oo 

= lim (f 0 (d° - a°) • • • / M (a M - a™)e iuAh f M (a M - af) • • • / 0 (a M - o")) eq | lt=0 (IV.30) 


Proof. By the Cauchy-Schwarz inequality and (IV.29), 

</o(a° - a2) • • • / A f(a M - a™)Ab 3 e iuh f M (d M - af) • • • / 0 (a M - af)) 
I (/o(a° - a°) ■ ■ ■ f M (a M - a**)Ab 2 f M (a M - a AI ) ■ ■ ■ f 0 (a M - a 7/ )) 


< 


eq 

1/2 


eq 

M _M\ 


(/o(a u - a2) • • • /m(S m - af )Ab*f M (a M - af ) • • • / 0 (a M - a")) 


1/2 


< an upper bound independent of N. 


This implies that the third derivative of the ‘characteristic function’ 

(/o(a° - a2) ■ ■ ■ /m(« M - a™)e zui f M (a M - a 7/ ) • • • / 0 (a M - a 77 )) eq (IV.31) 

is uniformly bounded on any bounded interval of u. Hence, by Theorem 4 and Lemma 6, Eq. (IV.27) is proved. □ 


As an example of a class of functions that satisfy the conditions of the above theorems, let us consider a Schwartz 
function which is defined as follows. 

Definition 1 (Schwartz function [8-11]). A Schwartz function is a infinitely differentiable function such that for any 
n,m = 0,1,..., 

sup\x n f {m \x)\ < oo. (IV.32) 

a?GR 


For example, x k e~ x is a Schwartz function for any k > 0. 

A Schwartz function has good properties, such as (i) a Schwartz function is absolutely integrable on R, (ii) the 
Fourier transform of a Schwartz function is also a Schwartz function. As a result, Schwartz functions satisfy all the 
conditions of the above theorems. 

Theorem 7. Assume that the conditions of Sec. IVA are satisfied for Aq, ..., Am,B. For any Schwartz function f, 
all the conditions of Theorems 5 and 6 are satisfied. 

Proof. By Lemma 4, there exist positive constants Kq, Ki and K 2 such that 


< 


</o(o° - a°) • • • f M (a M - a™)Ab 2 f M (a M - a. M ) • • • / 0 (a M - a. M )) eq| 

J d . . . dfi^M | /o[ £ ]/o [£/] . . . f M [Z]f M [f] | |( e *£o(o°-<*2) . .. e iZM(a™-a!<) Ah 2 ei P M (a™-a* 

Kq + K x £(|&| + $1) + K 2 J2 l»l 


, e i 5 '(a°-a 0 .)^ 


< / ^2^0 ... |/ 0 [3/ 0 [3] • • • /m[£]/m[£']| 


27r 


27T 


3,1 


Since i£jfj[£j] is the Fourier transform of /j(x), it is a Schwartz function. Hence, the integrand above is absolutely 
integrable, and Eq. (IV.26) is satisfied. Eq. (IV.29) is also proved similarly. □ 


Next, we prove theorems for the average value. 
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Theorem 8. If fj G L 2 (R) and f\fj [£] | 2 are absolutely integrable then 


lim [ da 0 , • • • da’f 1 (/ 0 (a° - a°) ■ ■ ■ f M (a M - a^)Abf M {a M - a • • • f 0 (a M - a^)) eq 

TV—>-oo J 

= [d §o • • • ^m|/oKo ]| 2 • •' I/mM 2 lim <e’« oAa ° • • • e ^A«“ A j e -ie M Aa M ... e -i£oAa°) (IV.33) 

J N—>oo 

Note that the conditions of this theorem are satisfied if, for example, fj are continuous and piecewise differentiable 
and fj (on each of their subdomains) are square integrable. 

Proof. 


J da 0 , ■■■ da^(f 0 (a° - a°)■■■ f M {a M - a^)Abf M (a M - a A J) • • • / 0 (a M - a^)) eq 

= J da, ■ ■ ■ da , 1 J ■ ■ ■ C ^^/o[£o]/o*&] • • • /mN/m[4] 

x ^ e *£o(Aa°-Aa°) _ _ _ e i$ M (Ao M -Ao,")^ e -i^(Aa"-Aof) . . . e -i^ 0 (Aa°-Aa°)^^ 

= J df o • • • d^lfoM 2 ■ ■ ■ |/M^M]| 2 (e^ Aa ° • • • e ^ A&M Abe-^ A&M ■ ■ ■ e~^ A& °) eq . 

By Lemma 3, there are 0(l)-positive constants K 0 and K\, and the integrand of the last integral is bounded by 


l/o[£o]| 2 • • • I/mM 2 (e i?oAa ° • • • e^ A&M Abe-^ A&M • • • e -^ Aa °) 


eq 


M 


< l/oKo]| 2 • • • |/m[£m ]| 2 ( Ko + I <1 ^2 l£j| 

1=0 


By assumption and the fact that the Fourier transform of any L 2 (R) function is also an L 2 (R) function, the above 
upper bound is absolutely integrable: 


J d£o ■ ■ ■ ^m|/o[^o ]| 2 • • • |/m[£m ]| 2 + K t |^|j <oo. 


Hence, by the dominated convergence theorem, 


lim / 

N—>oo J 

da° ■ ■ 

O • 

<3 

1 

o 

<C3 

^3 

) • • • fu(d M - a™)Abf M (a M - a. M ) • • • f 0 (a M 

a* 1 ))eq 

= lim 
N—fo o J 

[ dfo- 

'' df M |/o[£o] 2 ''' 

\fM^M]\ 2 {e iCoAa ° ■ ■ ■ e ^MAa M A l e -ii M Aa M . . . 

— i£oAa° \ 
c /eq 


■ ■ -d^M /o[£o]| 2 • • • /m[£ 

M \ I 2 lim t e iioAdP ■■■e iiMA&M Abe~ iiMA&M ■■■ 
N—foo 

-i£ 0 Aa°\ 
c /eq 


Theorem 9. If fj G L 2 (R) and /•,[£] are continuous and piecewise differentiable, and if (,fj[£\fj[£,\ and f 2 \fj[£,]\ 2 are 
absolutely integrable, then 

lim f da 0 --- daf Aa l .(f 0 (a° - a°) ■ ■ ■ f M {a M - a^)Abf M (a M - a,■ ■ ■ f 0 (a M - a^)) eq 

JV—>-oo J 

= jdf 0 ■■■ df M \foM 2 ■ ■ ■ f, .[ o - j ]! 2 irn(fimnmfi+i[^+i}\ 2 ■ ■ • i / mm 2 

x lim ( e i ?° Aa °... e i?«Aa M A j e -i?«Ai M . . . e -i{ 0 Aa°\ 

Ar->oo 

— ij d£o ■ ■ ■ cZ£m|/o[£o ]| 2 • • • |/m[£m ]| 2 

X lim — / e i ^ oAa ° • • • e *6-iAa !-1 e i(f 1 +s/2)Aa ! e i£ i+ iAa i+1 . . . g < Af Ao" 

N->oo ds ' 

X A be~ i ^ MA0,M ■ ■ ■ e -<i+ lA “ , + 1 e -i(6+s/ 2 ) Aai e —<(- lAa ' _1 . . . g-j£oAa°\ 


5—0 
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Note that the conditions of this theorem are satisfied if, for example, fj are continuously differentiable and piecewise 
twice differentiable, and (on each of their subdomains) are square integrable. 

Proof. 


J da 0 ,--- da ^ 1 A a l ,{f 0 (d° - a°) ■ ■ ■ f M (a M - a^)A bf M (a M - a* 1 ) ■ ■ ■ f 0 (a M - a^)) eq 

= J da 0 ,-■■ day A a[ j < ^° ■ ■ • ^^/o[£o]/o*&]''' /m[?m]/mKm] 

x ( e ?; £o(Aa°-Aa°) • • • e ^M( A a M - A a")A 5 e -*£M( A “ M - Aa “)... e -i £o( A “° -Aa *)) eq 

= i J ' Ul 'i ■ ■ ■ da-y j d -^° ■ ■ ■ ^ M 2 f M d{ ^ d _ ^ {/oN/o [&) ■ ■ ■ fM[tM]r M [?M] 

x (e^ Aa ° -Aa *) • • • e ^«( AaM — Ao “)ASe _ ^"( AaM_Aa * f ) • • • e _ ^^ Aa ° _Aa *)) eq | 

= -ij- s Jdf o • • • df M \f 0 M 2 ■ ■ ■ - |]i/ ;+ i[^+i]i 2 • • • i/m [£m]i 2 

x (e , ^° Aa ° • • • e i £ i - lZ ^“ i-1 e®(£ ! + s / 2 )'^“ , g i £ i + 1 '^“ !+1 ... 

x Abe~^ MA&M ■ ■ ■ e - i fi+ iA “ ,+1 e - i (i!- s / 2 ) Aai e - i fi- iAai_1 ... g _i ^ oAa °) eq | 

= j d £ 0 ■ ■ ■ d£M\fo[€o \\ 2 ■ ■ ■ |/;-iKi-i]| 2 I m (//[6]/;*Ki])l/+ifc+i]| 2 ''' I/m[£m]| 2 
x (g®£°^“ 0 ... g*?M^“ M /\^g _ '^M^“ M ... e~ i ^ oA ° , °) e q 

— i J d£o ■ ■ ■ (&;m|/o[£o]| 2 • • • |/m[£m]| 2 

x _/L/gi£oAa° _ _ _ e i{ i _ 1 Aa i_1 e i({ i +s/2)Ao i e iJ i + iAa i + 1 _ _ _ gif M Aa M 

ds ' 

x Abe~^ MAa,M ■ ■ ■ g _i ^+ lA “ i+1 g _ ff^ i + s / 2 )‘^“ i g _i ^- 1 ^“ i_1 ... g -i £oAa°\ 


By Lemmas 3 and 4, there are 0(l)-positive constants K 0 , K\ and iv 2 , and the integrand above is bounded by 

/ M 

IfoM 2 • ■ ■ JW.Ai i }\ 2 IM//&]/r&])l l/+i[ 6 +i ]| 2 • • • I/mKm ]| 2 Ko + K { Y \Sj 


j=0 


M 


M 


\foM 2 ■ • • |/m[Cm]| 2 K 0 + K x Y |^| + K 2 I Y & 


j=0 


\j =o 


By assumption, the above upper bound is absolutely integrable: 

[dfo ■ ■ ■ df M \foM 2 ■ ■ • iz-ii^i]i 2 \Mfiwnm i/+i[6+t]i 2 • • • i/mMi 2 I k 0 +k, y & 


M 


r I M I M 

+ df o • • • d£ M \f 0 M 2 ■ ■ ■ I/mNI 2 K 0 + K t Y ICl + d <2 Y 1 C 

J \ 3=0 \j =0 


3=0 


< 00 . 


Hence, by the dominated convergence theorem, the theorem is proved. 


□ 


V. DERIVATIONS OF THE MAIN RESULTS 


We now derive the main results of the paper using the above theorems. 
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A. Assumption on macroscopic systems and conditions on / 

We assume that conditions (C1)-(C4) of Sec. IV A are satisfied for all additive observables of interest, A 1 ,A 2 ,--- , Am, 
where M = 0(1). We assume this not only for equal-time operators A\(t), A 2 (t), • • • , Am( t), but also for operators 
at different times Ai(t\), A 2 (t 2 ), ■ • • , > where t\,t 2 ,--- ,tM denote the times of measurements, all of which 

are finite and independent of N. 

Let us discuss this assumption for the case where the system is composed of several units, each with a finite 
dimensional Hilbert space. For equal-time operators, the above assumption is believed to hold generally except in 
special regions, such as critical points, in the thermodynamic configuration space. For different-time operators, their 
correlations are bounded by the Lieb-Robinson bounds [12, 13] and therefore by the equal-time correlations, if the 
interaction is finite range. Hence, the above assumtion holds for many quantum systems except in special regions 
in the thermodynamic configuration space if the interaction is finite range and if the system is composed of several 
units, each with a finite dimensional Hilbert space. 

For the case where each unit has an infinite-dimensional Hilbert space, some of additive observables of interest might 
have infinite norms, for which condition (Cl) is not sastisfied. We have not yet succeeded in weakening condition 
(Cl) (while the Lieb-Robinson bounds have been extended to certain continuous quantum systems [14]). Therefore, 
at present, it is necessary to introduce a certain cutoff, by which each unit is well described by a finite-dimensional 
Hilbert space. If such an effective theory can be constructed conditions (C1)-(C4) can be satisfied. 

For the conditions on /, we will describe a sufficient condition on / for each result. Hence, a sufficient condition on 
/ for all the results is the logical conjunction of these conditions. We think the conditions are general enough from a 
physical viewpoint, although it might be possible to prove the results under weaker conditions on /. 

Note that Schwartz functions satisfy all the conditions on /, as is easily seen from their general properties and 
Thoerem 7. 


B. Derivation of Eq. (3) 

A sufficient condition on / for Eq. (3) is that / is a bounded piecewise-continuous real function. 

By the QCLT, the probability distribution of fluctuation of any additive quantity which satisfies the conditions in 
Sec. IV A converges weakly to the normal distribution in the TDL. Therefore, for any bounded piecewise-continuous 
real function /, 


p{a.) = lim ({/(a — a,)} 2 ) e 


N—> oo 


|/(a — «»)| 2 
( 27050 ^ 

\mf 


exp 


exp 


(a - (a)eq) 2 

2£< 

(x + A a.) 2 


26 a, 2 q 


J ( 277 ^) 1/2 

For example, in the case where / is gaussian f(x) = (2ttw 2 )~ 1 ^ 4 exp(— x 2 /4w 2 ), 


da 

dx. 


(V.l) 


P(a.) = 


^2w{w 2 + 6al q ) 


exp 


1 A al 


2 w 2 + Sa 2 q 


(V.2) 


C. Derivation of Eq. (4) and Eq. (13) 

A sufficient condition on / for Eqs. (4) and (13) is that / is a bounded piecewise-continuous real function and that 
conditions (IV.26) and (IV.29) are satisfied. 

To derive Eq. (4) and Eq. (13), we calculate the ‘characteristic function’ of A b(t) for the state p(a,): 


{exp(iuAb(t))) a , = —U-(/(a - a.) exp(mA6(t))/(a - a.)) 


P(a.) 


eq? 


(V.3) 
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where u is a real parameter. The expectation value and the variance of b(t) are obtained, respectively, by 

(Ht))a. = ~iS~ {exp(iuAb(t))) a 

OU I u=0 

2 o2 

- (kt))a.) )a. = ln(exp(mA6(t))) a . 


u =0 


Furthermore, by Theorems 5 and 6, we can interchange the limit and differentiation. Hence, by Theorem 4, 

{exp(iuAb(t))) a . 

1 1 


<W/[£]/[£' ]( e i f (a - a - ) e i “ AS(t) e I «'( a -°-)) eq . 


p(a .) 27t , 

1 1 


eq 


- (£ + ^)«(^{Aa, A6(t)}) eq - \u 2 5b(t) 2 eq 

-(€-0«(^[a,6(t)])eq] 


- 7 ~~r r—^ e u2sm ^' 2 [ dxf(x-u(^[a,b{t)]) eq )f{x + u(^[a,b{t)]) eci ) 
PM ^2v5a 2 eq J 2l 2l 


x exp — 


(x + Aa. — iu(^{Aa, A6(f)}) eq ) 2 ] 


2 Sal q 


^ pW) 6 U26Ht) ^ /2 y( a *~ iU (\{ A& ’ A k t )})eV U (^M0)])eq)- 

In the last line, we have defined the following function: 


If i 

Q{a'.,y)= / / dxf(x-y)f(x + y)ex p[--— -(x + Aa'J 2 

^2ir<y a 2 q J 2S < 


Note that g(a^,0) = p(a',) and {d/dy)q(a'„y )| y= o = 0. 
By using these relations, we have 

d 

(b{t)) a . = -* ln(exp(mA6(t))) 0 . 


n=0 


u—0 


d i l „ 

-—>-* w-lng(a. - m(-{Aa, A6(t)}) eq ,u( — [d,6(t)])eq) 

N—>oo OU Z Zl 

= -<^{Aa,A6(t)}) eq ^-lng(a.,0) - i(^[d, 6(t)]) eq hi q(a.,y) 
= -(^{Aa, Afe(<)}) eq (lnp)'. 


v=o 


Similarly, 


(m~m) a .f) a . 

d 2 

= - ln(exp(mA6(t))) a . 


•u=0 


<9 2 i i 

“ du? ln 9( - a * “ iu ( 2 { A “’ A6 ^^ eq ’2i & W]} eq ) 

,1 


= <H>; 


n—0 


{Aa, A6(t)})^ 


d 2 


eq \ 2 l-,“^J'eq^2 


lng(a.,0) 

5 a 


+ z(i{Aa,A6(t)})eq(^r[a,6(t)])eq ^-^lng(a.,y) 


j /=0 


- (^r[Aa, A&(t)]) 2 q ^lng(a.,y) 


y =o 


= ^eq + (^ { A a, A6(t)}) 2 q Inp(a .) - (^ [Aa, A6(f )]} 2 q Jh ln <?(a., j/) 


(V.4) 

(V.5) 


(V.6) 

(V.7) 

(V.8) 


(V.9) 


(V.10) 


y =o 
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D. Derivation of Eqs. (5) and (15) 


Since Eq. (5) can be obtained as a special case of Eq. (15), we prove the latter. 

A sufficient condition on fj for Eq. (15) is that fj £ L 2 (R) and £|/j[£]| 2 are absolutely integrable. This condition 
is satisfied if, for example, fj are continuous and piecewise differentiable and /j (on each of their subdomains) are 
square integrable. 

By Theorem 8, and /dx \fj{x )\ 2 = 1 and /dx x\fj(x )\ 2 = 0, and the spectrum decomposition of fj(a 3 — a 3 ,), 


A a 3 , = J da 0 ,da\ ■ ■ ■ da 3 , A a 3 ,(f 0 (d° - a°) • • • fj(a 3 - a 3 ,)fj(d 3 - a 3 ,) ■ ■ ■ f 0 (a M - a“)) eq (V.ll) 

= J daedal ■ ■ ■ da 3 ,~ x dx (/o(a° - a°) ■ ■ ■ /^(a- 7-1 - a 3 ~ 1 )f j {x) 

x (Ad 3 - x) fj(x)fj-i(a 3 ~ 1 - a 3 ,- 1 ) • • • f 0 {d M - a“)) eq (V.12) 

= J da°,dal ■ ■ ■ da 3 ,- 1 (/ 0 (a° - o^) ■ ■ ■ /j-i(d J_1 - a2 _ 1 )Aa J /.7-i(a J_1 - ■ • • fo(a M - a^)) eq 


-—> [ dfo • • • d£j-i\fo[£o \\ 2 • • • \fj-i[£j-i \\ 2 }im (e* A5 ° ... e ^- lA2J 'Aa^e 

iv—>-oo J N—>oo 


Aa J 


1 . . . g-*foAa°\ 


eq 


r 1_) 

= / o- • • dfj-i\f 0 [f 0 ]\ 2 ■ ■ ■ \fj-i[fj-i]\ 2 i y^g;([a ; ,a J ])e q 

1=0 


= 0. 


(V.13) 


E. Derivation of Eq. (7) and Eq. (16) 


A sufficient condition on / for Eqs. (7) and (16) is that fj G L 2 (R) and fj[f] are continuous and piecewise 
differentiable, and that £/j [£]/)[£] and f 2 \fj[f ]\ 2 are absolutely integrable. These conditions are satisfied if, for 
example, fj are continuously differentiable and piecewise twice differentiable, and fj , xfj , /j, f" (on each of their 
subdomains) are square integrable. 

We first prove Eq. (16). By Theorem 9, we have for j < k 


Aa{Aa^ 

= J da°,dal • • • da k , Aa 3 ,Aa k (f 0 {d° - a°) • • • f k (a k - a k )f k (d k - a k ) ■ ■ ■ f 0 {d M - a")) eq (V.14) 

= J da°,dal ■ ■ ■ da k ~ x Aa 3 ,(f 0 (a° - a°) • • • f k - i(a fc_1 - a k ~ 1 )Ad k f k -i(a k - 1 - a^ -1 ) • • • f 0 (a M - a^)) eq 

J df o • • • d£fc_i|/o[£o]| 2 • • • l/j-i[6-i]| 2lm (/j[Cf]/j Kf])l/j+i[Ci+i]| 2 ■'' \fk-i[fk-i\\ 2 

x lim (e^° Aa ° • • • Ad fc e _ ^ fc_lA “ fc_1 • • • e -*°Aa<\ 

N—too 

+ j d£o ■ ■ ■ d6c-i|/o[£o]| 2 • • • |/fc-i[£fc-i]| 2 

X lim Re(e^° Aa ° • • • ( 'A A "' A " AdVA • -W ... ^ k ^Aa k ~ 1 

JV-> oo ' 


x Ad k e~ iik ~ lA&k 1 • • • e -*?° Aa °) eq _ (V.15) 
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By Theorem 3, we can evaluate these limits by the QCLT as 

J d £o • • • d£fc-i|/o[£o]| 2 • • • \fj-i[€j-i ]\ 2 Im (fj ] // [£j ]) | /j+1 [£j+1 ] | 2 • • • \fk-i[€k-i ]\ 2 

lim (e i|oAa ° • • • • • • e~ j;&Aa ° 

N—> oo 


} eq 


• • • d^\f 0 M 2 ■ ■ •i/ J -ife_i]i 2 im(/'fe]/;[o])i/ J+ ife + i]i 2 • • ■\f k ^k-i]\ 2 j2^i(-[d l ,d k ]y 

i=o 1 

= ~ j [0])^(v[a i ,a fe ]) eq 

= (v[a- 7 ,a fc ]) eq Im J dxxf j (x)f j (x) 

= 0, 



! di o • •• ^ fc -i|/oKo]| 2 • • • |/fe-i[6-i]| 2 lim Re(e i « oAa ° • • • , :A " V- A " AuV- w " ■ ■ ■ e^- lAa ' 

J N=r oo 

Aa fc e _i?fc_lAa&_1 • • • e" ifoAa °) 

( fc—1 j—1 \ 

{cP, a fc }) eq + ^<7^*. « fc ])eq E aJ ']>e 

= (1 {a^a fe }) eq + ^(i[a',a fc ]) e 2 q |d& & 2 |/,&]| 2 


;=o 


h =0 


1=0 

1-1 


= (J{a^a fe }) eq + ^F / (^[a^a i ]) eq (i[a i ,a fe ]) 


;=o 


2i 


Here, 

F ; = 4 J dx{fl(x)} 2 , 


where f{(x) denotes the Fourier transform of *£/;[£]■ For the case where fi(x) is twice differentiable, Fi 
rewritten as Fi = —4 f dxf”(x)fi(x). 

We finally prove Eq. (7) using Eq. (16). Putting a 0 = a and a 1 = 6, we have 

Ac^Aai = J da°,dal Aa^Aoi(/ 0 (a° - i) - aDf^aftti) - ai)/ 0 (a° - a°)) eq 

= J da°,dx Aa“{/i(i)} 2 (/ 0 (a° - a°)(Aa 1 (ti) - x)f 0 (d° - a°)) eq 

= J da 0 Aa°(/ 0 (a° - a°)Aa 1 (ti)/ 0 (a° - a°)) eq 
= Ao 2 (Aa 1 (t 1 )) a o. 


VI. RESULTS FOR THE CASE OF GAUSSIAN f(x) 


A typical f(x) is gaussian, 


fi x ) = , 0 L/4 ex P - 

(27TUA) 1 / 4 

with w = 0(1) > 0. In this section, we present results for this case. 
The measurement error is simply given by 


x 2 
4 w 2 


(V.16) 


(V.17) 


(V.18) 
can be 


(V.19) 


(VL1) 


(VIA) 


6a e rr = W. 
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Eq. (4) reads 


Eq. (13) reads 


( A Ht))a, 


m 


(|{Aa, Afr(f)}) eq 
(w 2 + Sa^ q ) 


(VI.3) 


m)-m) a .) 2 ) a .-5b 2 e = 


w 


-1 

2 -L Sal 


[UAa,Ab(t)}) 


eq 


W 


2 '2i I 


*> &(*)]> 


2 

eq* 


In Eq. (16), and Fj are given by 


(VI.4) 


for 


Sa J er 2 r = tu|, (VI.5) 

^ = V^ 2 (VI.6) 


/jO) 


(27rw|)i/4 6XP 


X 2 


(VI-7) 
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